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Abstract 

We consider a classically scale invariant version of the Standard Model, extended by an extra 
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1. Introduction 


In 2012, the ATLAS and CMS experiments of the Large Hadron Collider (LHC) at CERN reported 
the discovery [1,2] of a boson that appears to be [3,4] the long-sought Higgs particle [5-8] of the 
Standard Model (SM) [9-11]. The latest precise measurements of the Higgs boson mass [12] place its 
value at Mh = 125.09 ± 0.24 GeV . With the last piece of the puzzle now in place the SM seems 
complete. 

Although the SM is a perfectly consistent quantum held theory and is expected to be valid up to 
energies of the order of the Planck scale ( Mp ), where quantum gravitational corrections are believed 
to come into play, with the observed value of the Higgs boson mass the Higgs quartic self-coupling 
runs to negative values at an energy scale lower than Mp. This is the SM vacuum stability problem. 
Furthermore, if one considers the SM as an effective held theory embedded in a more fundamental 
theory involving other scales, then the issue of the smallness of the Higgs mass in comparison to Mp 
or other scales (like Mq xjt) arises and, as a result, the so-called hierarchy problem. In a general UV 
completion of the SM the Higgs mass would not be protected against contributions from the new 
massive states and fine-tuning would be needed. If the SM is embedded in a supersymmetric theory, 
the Higgs mass is radiatively stable down to the scale of supersymmetry breaking, which has to be 
in the vicinity of the electroweak breaking. Nevertheless, the new degrees of freedom would have to 
be observed in this neighborhood and the results from LHC have been negative for supersymmetry 
so far. 

Another possibility is that of scale invariance. In the case that at very high energies the Higgs 
mass vanishes and the theory, having no dimensionful parameter, is classically scale invariant, no such 
mass term can arise by radiative corrections at lower energy scales [13]. Only logarithms multiplying 
the tree-level terms can arise and the Higgs effective potential will have the form 

Ves(H) = A |IL| 4 + C\H\ 4 ln(|iL| 2 / p 2 ). (1.1) 

However, a dimensionful parameter and, in consequence, physical masses can arise through spon¬ 
taneous symmetry breaking. The associated emergence of a dimensionful parameter from a scale 
invariant theory has been termed dimensional transmutation and its realization through radiative 
corrections Coleman-Weinberg mechanism [14]. Minimization of (1.1) gives (H) ~ fie~ x ^/ c . Tak¬ 
ing the renormalization scale /i at the UV, we can in principle account for an exponential hierarchy 
between the electroweak scale ( H) and the Planck scale /r = Mp [15-18], provided the particle con¬ 
tent of the theory is such that the radiative corrections coefficient C is positive. Nevertheless, in the 
minimal version of the Standard Model the dominant contribution in C comes from the top quark 
and makes it negative, therefore no hierarchy is generated. This can change if extra bosonic degrees 
of freedom with sizable couplings are introduced. These could render C positive and sufficiently large 
in order to generate the desired hierarchy. The extra degrees of freedom could be scalars or gauge 
bosons, with the latter possibility corresponding to an enlargement of the gauge group. 

The observed Higgs mass Mh = 125.09 GeV, resulting in A(M*) « 0.1285 [19,20], leads to negative 
values for the Higgs self-coupling A(/r) above scales of 0(1O 10 GeV) which means the vacuum is 
actually metastable [21-31]. This behavior can be avoided in extended versions of the Standard Model, 
the simplest of which consists in the introduction of a gauge singlet scalar field that couples only to 
the Higgs field. One way that the extra field affects stability is through the positive contribution to 
the renormalization group equation of A induced by its portal coupling to the Higgs that tends to 
counterbalance the negative top quark contribution. Barring the unappealing possibility of introducing 
an additional ad hoc mass scale associated with this singlet, one may incorporate the new held in 
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the classically scale invariant version of the Standard Model [32-38]. This allows for a second way to 
affect stability due to the increase of the low-energy value of the Higgs coupling caused by the extra 
contributions from the vacuum expectation value (vev) of the singlet field. 

The new scalar can acquire its vev through dimensional transmutation by a realization of the 
Coleman-Weinberg mechanism (CW). Loop corrections to its effective potential can generate for it a 
nonzero vev that is fed to the standard Higgs through the portal coupling. However, as we remarked 
above discussing (1.1), the coefficient C of the one-loop effective scalar potential has to be positive. 
This can be readily achieved if the new field is charged under a new gauge group, since the new gauge 
boson contributions to C are positive. Thus, introducing an extra dark gauge group [39-55] with 
extra scalar degrees of freedom can in principle lead to a successful realization of CW. 

The new sector associated with UV stabilization and low-energy symmetry breaking can also in 
principle account for dark matter. A characteristic example is provided by dark matter being a 
vector boson, resulting from the breaking of an extra SU(2)x gauge symmetry [47,48], with mass 
in the TeV range. Depending on the model, additional dark degrees of freedom may be needed in 
order to saturate the dark matter relic abundance. These could be additional dark scalars [56-59], 
although fermions are allowed too if their masses are light enough not to invalidate the dark sector 
CW mechanism [60-62]. 

Regarding neutrino masses, we can introduce right-handed neutrinos as singlets [63-70], neutral 
both under the Standard Model as well as the dark gauge group. A mass term for them can arise 
from a coupling to a total singlet scalar that obtains a vev through its couplings to the dark sector. 
Thus, neutrinos may acquire masses as a result of a seesaw mechanism driven by the scale of the dark 
sector. Nevertheless, this low-energy seesaw procedure still requires small Yukawa couplings for the 
left-handed neutrinos of the order of the corresponding electron Yukawa coupling since the singlet 
vev, tied to the gauge symmetry breaking, cannot be too large. 

In the present paper we have reconsidered an SU(2)x extension of the classically scale invariant 
(CSI) version of the Standard model and, in the light of the breaking of the gauge symmetry SU (2) l x 
U(l)y x SU(2)x —> U( l)em, we address the issues of UV stability, the structure of the multi-Higgs 
sector, dark matter and neutrino masses. Apart from the right-handed neutrinos and a real scalar 
singlet, having a coupling to them and to the standard Higgs, the model consists of a dark SU(2)x 
sector composed of three dark gauge vectors and a dark isodoublet possessing portal couplings to 
the standard Higgs and the aforementioned singlet. After symmetry breaking, the extended Higgs 
sector consists of three states, one to be identified with the standard Higgs and two additional scalars. 
The resulting dark vectors, being stable and weakly interacting massive particles (WIMPs), can be 
identified with dark matter. For a suitable range of the free parameters of the model all relevant 
experimental constraints can be met as well as stability and perturbativity. In addition, the model 
has a range of definite dark matter predictions which can be either falsified or verified by new limits 
or observations in near future experiments. 

The paper is organized as follows. In the next section, we present the model and analyze the 
stability of the tree-level potential. We proceed to obtain possible flat directions, setting up the 
model for the study of symmetry breaking through the Coleman-Weinberg mechanism. We compute 
the one-loop effective potential and the resulting scalar masses. Subsequently, in Sec. 3, we undertake 
a phenomenological analysis of the model. We identify one of the predicted scalar states with the 
observed Higgs boson. Then, we find benchmark sets of values for a minimal subset of the free 
parameters of the model that correctly reproduce the Higgs boson mass. After that, we scan over the 
rest of the parameters and obtain masses for the dark gauge bosons, the right-handed neutrinos and 
one of the scalar bosons, all the while checking that the stability and perturbativity constraints are 
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satisfied. In Sec. 4, for the same set of benchmark values, we calculate the dark matter relic density 
and constrain the masses of the dark gauge bosons from both the observed relic density and the limits 
set by direct detection experiments. Finally, in Sec. 5 we summarize and conclude. 


2. The model 

In this section we present the model and study its properties. Employing the Coleman-Weinberg 
mechanism [14] in the Gildener-Weinberg [71] formalism we minimize the tree-level potential and find 
the flat direction between the vevs of the scalar fields. Then we obtain the tree-level masses of the 
scalars, one of which (it may be called darkon) turns out to be massless due to the flat direction. 
Including the one-loop potential, we find that radiative corrections become dominant along the flat 
direction and lift the darkon’s mass to values that can be even higher than the masses of the other 
scalars. 


2.1. The tree-level scalar potential 


In order to address the open issues discussed in the Introduction we consider the Standard Model in 
a classically scale invariant (CSI) framework and extend the gauge group with an additional SU{2)x 
symmetry [47,48]. In addition to the new gauge bosons, the dark sector contains a SM-singlet scalar 
SJJ{2)x isodoublet 4>. Aiming at the problem of neutrino mass generation, we also introduce a real 
scalar a, singlet under both the SM and dark gauge groups. Right-handed neutrinos are also included 
in the standard fashion as total fermionic singlets. The tree-level scalar potential, in terms of the SM 
Higgs field H and the new scalars <!>, a, has the form 

Ho = A h (H^H) 2 + A^^) 2 + - A hcj> (H^H){^^) - + ^(H^H)a 2 , (2.1) 

where we included all possible couplings among scalars and have introduced negative signs for the 
portal couplings A h<f> and \<f, a . Models having all mixing scalar couplings positive, i.e. both the portals 
to the dark sector A h<f>, A ^ and the observable sector Higgs mixing A ho, do not lead to a symmetry 
breaking flat direction and are not of interest. It is then reasonable to examine models with the portal 
to the dark sector being negative. Taking the observable mixing A ha also negative is not necessarily 
interesting because it allows for a flat direction independent of the dark sector. Therefore, we restrict 
the possible breaking patterns making the above choice of signs which is sufficient for our purposes. 

In addition to the scalar potential, the Lagrangian has the following extra Yukawa terms: 


- C N = Y^li ia 2 H*Nj + H.c. + Y^NfNja, (2.2) 

where is the Dirac neutrino Yukawa matrix which couples the left-handed lepton doublet L; to 
the SM Higgs doublet H and the right-handed neutrino Nj and Y& is the right-handed Majorana 
neutrino Yukawa matrix which will be assumed diagonal. 

Considering the unitary gauge and the symmetry breaking pattern SU{2)l x U{1)y x SU{2)x —>■ 
SU{2)l x U{1)y U{ l)em, we may replace the scalar doublets by 






Then, the tree-level potential takes the form 

V 0 (h.i,,a)=L h * + L^ + L^ 


YftV _ + YfcV 2 . 

4 4 4 


(2.3) 


(2.4) 
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The scalar potential is bounded from below if the matrix 
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2A h 

^h(f) 
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(2.5) 
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2Acr 

/ 



is copositive , i.e. such that ri a A a bVb is positive for non-negative vectors in the basis (h 2 , cj) 2 , a 2 ). It 
can be shown [72-76] that this is equivalent to the conditions 
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(2.7) 

A/lCT 


2 A^A (j 21^/ 


Note that the last condition is equivalent to either of the following statements: 


+ 


V<r i 

A—V — 


2-^/A^A^ 2\/A7A7 2^/A^Ao- 

det *4, — AhA^A ct 4 ^A^Act T A^A^, A ( , )(T A/ ; T jAhj/AherAj/cr — O' 


2-v/AhAcr 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


Therefore, vacuum stability requires the validity of the above conditions to hold at all energies up to 
Mp. In order to study the flat directions of the tree-level potential we may parametrize the scalar 
fields as 

h = (pNi, 4> = tpN 2 , (J = (2-11) 

with iVj a unit vector in the three-dimensional held space. Then, the tree-level potential attains the 
form 

4 

Vb = ^ [AhiVi + A^JVf + ActM 3 4 " Ah^iiVl + A^^i N 2 - \^N 2 N 2 ] . (2.12) 

The condition for an extremum along a particular direction N t = m is [71] 


5F 0 

dNi 

= Vo(n) = 0. 

n 

(2.13) 

Then, the equations giving the symmetry breaking direction are 


2A h n\ 

— ^h(j) n 2 ~ A her ^3 

(2.14) 

2A^77-2 

— A/j^TT--^ + A^^Tlg 

(2.15) 

2A CT n| 


(2.16) 

A h'lT’l + A077.2 + Ao-Ug — 

1^2 A ( ^ cr 77-2^'3 + A/ icr ?'i-L??-g — 0. 

(2.17) 
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The solution of these equations in terms of the scalar couplings is 


n\ 


nk = 


rio = 


_4A tr Aj, A^_ 

^Afj^^A^-t- A A^^ ) A ha (2 A^-TA^fj ) 

2Acr (2A(£ + \(ptj ( A/j,^> A^>cr ) ^hcr (2A^+ A^ 

^h(f> ^cfxr 2\<p\hcr 

2Acr (2A0 + \h(())+ ^(f)(T ( ^h<f> A ) A her 2 A^> ~J - A (j)(j ) 


(2.18) 


Note that n\ + n\ + n| = 1. 


2.2. The scalar masses 

Assuming spontaneous breaking of the gauge and the scale symmetry via the Coleman-Weinberg 
mechanism, we can write the shifted scalar fields as 


h = (ip + v)ni, 4 > = (</? + v) n 2 , a = (ip+ v)n 3. 

The individual vevs are 

(h) = v h = vn 1, (</>) = v ( j ) = vn 2 , (a) = v a = vn 3 . 

From the shifted tree-level potential we can read off the scalar mass matrix 

2A h n\ -nin 2 \ h(j> +mn 3 \ ho 
Mq = v z ( -?rin 2 A/ l 0 2A^n2 -n 2 n 3 A0 o 
+nin 3 X ha -n 2 n 3 \q, a 2 \ a n 3 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


in the (h , q i>, c) basis. We can now set up the diagonalization of the mass matrix (2.21) by introducing 
a general rotation in terms of three parametric angles, 


nMln -1 = Ml 


with the rotation matrix 7Z 1 given by 


cos a cos /3 


sm a 


cos a sin f3 


1Z = I — cos f3 cos 7 sin a + sin (3 sin 7 cos a cos 7 — cos 7 sin a sin /3 — cos /3 sin 7 
— cos 7 sin (3 — cos j3 sin a sin 7 cos a sin 7 cos (3 cos 7 — sin a sin (3 sin 7 


and 


= . 7Z- 


o 



( 2 . 22 ) 


(2.23) 


(2.24) 


Next, we choose two of the above three angles in the rotation matrix to parametrize the total vev v 
direction according to 

Vh = v sin a = vn\ 

V(j, = v cos a cos 7 = vn 2 (2.25) 

v a = v cos a sin 7 = vn 3 . 
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Then, M. 2 d is diagonal, provided that the following relations are satisfied: 


9 vt 

tcLUA Qj - H/ _ _ _ 

V c()~^~^cr 2^ A ct Xh<f> A <£Xhcr ) ^h<r) 


4A^>A tT -A| cr 


tan 2 7 = 4 = 2A ^ 0 T \t Xha 

v a. 4Afr Act — X ha 


(2.26) 


t _ T , n,q _ »fet7ty »(A hCT +A h ^,) 

(A^+A ct +A^ ct )^);2 -A, 1 ^d 2 ' 

The resulting mass eigenvalues are 

2 / 222 2 / \2 

M d j2 =\ h vj l cos a cos /3 + X^v^ (cos (3 cos 7 sin a — sin (3 sin 7 ) 

+A ct u 2 (cos 7 sin /3 + cos f3 sin a sin y ) 2 

+Xh cos a cos (3 (cos (3 cos 7 sin a — sin (3 sin 7 ) (2.27) 

—X^v^Va (cos (3 cos 7 sin a — sin / 3 sin 7 ) (cos 7 sin f3 + cos (3 sin a sin 7 ) 

—Xh a VhV a cos a cos /3 (cos 7 sin (3 + cos (3 sin a sin 7 ) 

M 2 2 = 0 (2.28) 

2 / 222 2 / \2 

M ^ 3 /2 = A hVh cos a sin f3 + X^v^ (sin (3 cos 7 sin a + cos /3 sin 7 )" 

+X a v 2 (cos 7 cos (3 — sin /3 sin a sin 7) 2 

+A hxpVhVif) cos a sin (3 (sin [3 cos 7 sin a + cos f3 sin 7 ) (2.29) 

+A 'fxyV'fyVa (sin /3 cos 7 sin a + cos /3 sin 7 ) (cos 7 cos j3 — sin (3 sin a sin 7 ) 

+Xh a VhV ( j cos a sin j3 (cos 7 cos (3 — sin /3 sin a sin 7 ) 

As expected, one of these masses turns out to be zero at tree level. Regarding the rest, Mf l] and 
M hs are ultimately functions of the overall vev v and the scalar couplings. Since their exact, analytic 
expressions are not necessary, we leave them as they stand. 


2.3. Neutrinos 


One of the defining properties of the present model is that it incorporates the appropriate structure 
for massive Majorana neutrinos. Right-handed neutrinos in three families are introduced as singlets 
of both the Standard Model and the SU(2)x dark sector. They obtain their mass as a result of broken 
scale invariance through their coupling to the singlet cr that mediates between the two sectors and 
obtains a nonzero vev. The Yukawa terms (2.2) that give rise to neutrino masses are 


yb' 

-^=v h Vi ia 2 Nj + H.c. + Y?v a N i c N j . (2.30) 


The neutrino mass matrix, being of the seesaw type, can lead to the desired scale of 0(0.1 eV) for 
the left-handed neutrino masses. Thus, in a (y s , Nf) basis, we have 


0 

\ 72 Vh ) 


(2.31) 
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Assuming v\ l to be no more than the lightest charged lepton mass, namely 0(1O~ 1 * * 4 GeV), and 
taking characteristic values v a ~ 0(1 TeV) and Y a ~ 0(0.1), we have Y u Vh -C Y a v a and we arrive 
at approximate eigenvalues 


M n Rj Y* j V c r, m u 


2 

JY_y(ik) ( v -l\(kf) v {tj) 
\v a v 4 ’ v ’ 


(2.32) 


with M tv ~ 0(100 GeV) and m v ~ 0(0.1 eV). As we will see next, the right-handed neutrino 
mass scale is related to the masses of the rest of the particles and cannot take arbitrary values. 


2.4. The one-loop potential 

Now, let us consider the full one-loop potential. Following the Gildener-Weinberg approach [71] to 
symmetry breaking we have 
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wN + v ^ 


= 0 


#=u(n+<5n) 


or, expanding and using (2.13), 
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1 dNidNj 
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The last statement is equivalent to 


dV l 
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+ n-i aAT 

N=n 

N=n 

n% dNi 
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N=n 


9V i(ny?) 


d<p 
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(2.33) 


(2.34) 


(2.35) 


(2.36) 


In this approach the couplings of the tree-level potential are assumed to depend on the renormalization 
scale p and the tree-level minimization condition to be realized at a particular value p = A. Thus, 
along the minimum flat direction at the scale A the one-loop effective potential has the form 


Vi (nt£>) = Aip 4 + Bp 4 log ^ • 


(2.37) 


The coefficients A and B are dimensionless parameters and are given in the MS scheme by 
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(2.39) 


B = 


1 


647r 2 u 4 



Af£ + 6 M^ + 3M| + 9- 12Af t 4 



The condition (2.36) gives 



Thus the one-loop effective potential becomes 


Vi(n <p) = Bip A 



1 

2 


(2.40) 


(2.41) 


The pseudo-Goldstone boson (darkon) mass is now shifted from zero to 


2 _ d 2 Vi{nip) 
Mh2 “ d<p* 


(p=V 


^ (A7 4 + Ml + 6 Ml + 3 M% + 9M 4 - 12Af 4 - 6A7 4 ) , (2.42) 


where we have assumed for simplicity that all right-handed neutrinos are degenerate in mass. Note 
that the right-handed neutrino contribution to (2.42), as fermionic, enters with a minus sign. As a 
result, Afjv cannot be too large. In what follows we shall identify the state h\ with the observed Higgs 
of 125.09 GeV and choose a higher value for the /13 mass. The /12 state, although massless at tree 
level, can have any mass with respect to the other scalars due to the sizable one-loop correction. Note 
that radiative corrections to the tree-level masses of h\ and /13 are small enough to ignore to a first 
approximation. 


3. Phenomenological analysis 

In this section we present an analysis of the phenomenological viability of the model, taking into 
account theoretical and experimental constraints. Our procedure in broad terms will be as follows: 
First, we choose values from a subset of the free parameters of the model (i.e. v^, v a and some of 
the scalar couplings), appropriate to fix the mass A#/,, to the experimental value of 125.09 GeV. 
Then the Mh 3 mass is automatically obtained. In order to calculate the darkon mass Mh 2 we scan 
over the two remaining unknown masses Mx and Mjy in (2.42), while checking that the stability and 
perturbativity conditions are satisfied. Finally, we calculate the total decay rates of all the scalar 
bosons and compare the one corresponding to the Higgs boson with the bounds set by LHC. 

3.1. Theoretical constraints 

The tree-level potential (2.4) and the one-loop effective potential (2.41) have to be bounded from 
below for the vacuum to be stable. For this to be valid, the stability conditions (2.6)-(2.10) need to 
hold for all energies up to the Planck scale (Afp = 1.22 x 10 19 GeV) as well as the positivity condition 
B > 0 has to be satisfied. The latter translates to 

Ml + 9 Ml - 6 Afy > 12Af 4 - 6 M^ - 3Aff - (3.1) 

or 

Ml + 9 Ml - 6A4 > (317.26 GeV ) 4 , (3.2) 

where we used the values M t = 173.34 GeV [77], M\y = 80.384 GeV and Mz = 91.1876 GeV. The 
above inequality implies that the masses of the extra gauge bosons Mx have to be in general larger 
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than the masses of the right-handed neutrinos Mjv, unless the scalar boson mass Mj l3 is considerably 
larger than the right-hand side of (3.2). 

Another constraint arises from the requirement that the model must remain perturbative all the 
way up to Mp. This can be achieved by demanding that all couplings are bounded, 

all couplings <27r. (3-3) 


To determine how the couplings of the model vary with energy, we need to solve the renormalization 
group equations (RGEs). We present the two-loop gauge and one-loop Yukawa and scalar RGEs below 
(however in our numerical analysis we use the full two-loop RGEs for all the couplings, computed 
using Refs. [78-82]): 


Pgi = ^9? + ^?( 199 S? + 135 <72 + 440^ - 85y 2 ) 

Pgi = + ^2 <^2 ( 27 Si + 17532 + 3603! “ ^vf) 

Pgs = -7.9.3 + ^2^3.3(11.9? + 453 2 2 - 26O3! - 203?) 

„ 43 o 1 259 . 

P 9 x -~y 9x ~ {^y~^ 9x 

Pyt = m ( 2 y* ~ 20 91 ~ I 92 ~ 893 ) 

Pv a =AY a Tr(Y a Y*^ + 12Y a Y*Y a 

Px h = ~^Vt + 24A! + A h ^12 y 2 - -g\ - 9 g^j + i^9i + ^ 9 i 92 + + 2A+ -A \ a 

Px,4, = g 9 x “ 9 dx^<f> + 24 A| + 2A)j0 + - X^ a 

Px„ = - 64Tr ( Y„ Y* Y„ Y* ^ + lOA^T^Y;) + 18A* + 2\ 2 ha + 2X 2 a 

( 9 9 9 9 9 9 

Px htj , = A h<p ( Gi/t + 12A/, + 12A <f, — 4A^0 — — g 1 — -g 2 — -3x ) + 'W'V< T 

Px^ = \t>a (^8Tr (Y a Y*^j + 12A^ + 6A CT — 4A0 CT — — g\^j + 4A/ lcr A/i0 

Px ha = ^ha (ftyi + 8 Tr(h (r Y*^j + 12 A/j + 6 A<j + 4A/j CT — —y 2 ~ + 4A/ l< ^ ) A^ CT , 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


where we defined /3 K = ( 47 r ) 2 3 ^ 7 . 

In order to solve the RGEs we have to specify the boundary conditions for the couplings. For the 
SM gauge couplings and the top quark Yukawa coupling we use the NNLO values at M t [22,83]: 
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(3.16) 


5i (M = M t ) 
.92(5 = M t ) 

93 (/i = M t ) 

5t(+ = M t ) 


- ( 0.35830 + 0.00011 


- 173.34 - 0.00020 


0.64779 + 0.00004 


M t 

GeV 


Mi 

GeV 


- 173.34 +0.00011 


f M w 


- 80.384 GeV 




1.1666 + 0.00314 (_ o 00046 
V 0.0007 J 


0.014 GeV 
AIw — 80.384 GeV\ 
0.014 GeV J 

M t 


- 173.34 


0.93690 + 0.00556 


M t 

GeV 


- 173.34 - 0.00042 


GeV 

a s (M z ) — 0.1184A 


0.0007 




(3.17) 

(3.18) 

(3.19) 


In our numerical analysis we use as inputs the central values a s (M z ) = 0.1184, M\y = 80.384 GeV 
and Mt = 173.34 GeV. For the right-handed neutrino Yukawa coupling Y a and dark gauge coupling 
gx we define Y a {Mx) = Mx/v a and gx(Mx) = 2 Mx /'(+ respectively. Finally, to define the scalar 
couplings A i, we consider their values at the renormalization scale A determined by (2.40), where the 
one-loop effective potential is minimized [84]. 

A few comments are in order regarding the behavior of the running couplings. First of all, the 
SU(2)x gauge coupling gx decreases at higher energies, being asymptotically free, in a similar fashion 
to the SU(2)l gauge coupling g 2 . On the other hand, the RGE of the right-handed neutrino Yukawa 
coupling Y a has a positive sign and forces Y a to increase with energy until it potentially reaches a 
Landau pole. As it turns out, this can be avoided if Y a {Mx) < 0.35. The Higgs self-coupling A/, 
generally behaves like the corresponding one in the SM. There, Ah drops fast at increasing energy due 
to the large negative contribution from the top Yukawa coupling yt , crosses zero at some point and 
then becomes nearly constant up to the Planck scale. Nevertheless, in our case, we have the freedom 
to choose a starting value for A/, such that it remains positive inside the whole energy range under 
consideration. The self-coupling of the singlet scalar A a depends highly on Y a and it too can reach 
a Landau pole unless Y a (Mx) < 0.31. Therefore Y a is further constrained. Now, the dark scalar 
self-coupling A^ generally increases with energy, driven mainly by the first term in (3.11). A Landau 
pole is avoided if we have gx(Mx) < 2.51. Finally, the scalar portal couplings (A/, 0 , A^, A/, CT ) are 
mainly multiplicatively renormalized and, as it turns out, they do not run much if we choose initial 
values that are small enough. 

The model contains many free parameters, therefore we need to restrict or fix most of them. The 
initial eight dimensionless free parameters 2 Ah, A^, A a , A/,. 0 , A ha, A 0 CT , gx, Y a are reduced to six after 
imposing the experimental values on Vh and M/,,. Whether we use the dimensionless scalar couplings 
or, alternatively, the vevs through the minimization conditions (2.14)-(2.16) as input parameters, is a 
matter of choice. A six-dimensional parameter space is not easily managed in its full generality. So, we 
propose to proceed in the following way: We leave gx and Y a free and list characteristic (benchmark) 
values for the scalar couplings (at A) and the vevs v a and that reproduce the measured Higgs mass 
Mh x = 125.09 GeV. These are shown in Table 1 where we also show the value for the mass of h .3 
that we obtain. 

2 We have assumed that the neutrino Dirac Yukawa coupling takes up values in the neighborhood of the corresponding 
electron Yukawa coupling. 
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Set 

Vh[ GeV] 

v$[ GeV] 

v a [ GeV] 

A/,, (A) 

A</>(A) 

A ct (A) 

Aft0(A) 

A0o-(A) 

Aft CT (A) 

M hs [ GeV] 

A 

246 

2112 

770 

0.1276 

0.004 

0.2257 

0.0036 

0.06 

0.001 

550.62 

B 

246 

3245 

1470 

0.1285 

0.0005 

0.0122 

0.0015 

0.005 

0.0001 

251.93 

C 

246 

4513 

2181 

0.1287 

0.0035 

0.0642 

0.001 

0.03 

0.001 

868.15 


Table 1: Benchmark sets of values for the model parameters able to reproduce the observed Higgs 
boson mass Mf l] = 125.09 GeV. 


Not all these sets are compatible with the stability of the potential. For example, with the values 
in set B the stability condition (2.10) is violated. Choosing the first set of values (A), we scan over 
the remaining two free parameters gx and Y a and obtain the darkon mass Mft 2 contours shown in Fig. 
1 , while checking that the stability conditions ( 2 . 6 )-( 2 . 10 ) are satisfied, that the one-loop potential is 
bounded from below (3.2) and that all the couplings remain perturbative up to the Planck scale (3.3). 


250 



500 1000 1500 2000 2500 


M x [GeV] 

Figure 1: (color online). Parameter space scan in the plane (gx, Y a ), taking into account constraints 
from stability and perturbativity. The color coding signifies the mass of the darkon Mft 2 ■ 

We also present the running of the scalar couplings in Fig. 2, again for the values of set A in Table 
1 and indicative values for gx and Y a corresponding to Mx = 725 GeV and Mx = 240 GeV. 
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M n = 240 GeV, M x = 725 GeV 



/i [GeV] 


Figure 2: (color online). The RG evolution of the scalar couplings at two-loop order for M x = 
240 GeV and M x = 725 GeV. 


3.2. Experimental constraints 


The three scalar fields of the present model all develop a vev and mix through the portal terms in the 
scalar potential (2.4). Moreover, the corresponding mass eigenstates interact with the SM electroweak 
sector, as well as with the SU(2)x gauge fields and the right-handed neutrinos. The strength of these 
interactions is suppressed by the corresponding entries in the rotation matrix 7Z (2.23). In order to 
study possible signatures of the extra scalars at the LHC and future colliders, we may construct an 
effective Lagrangian that contains all the interactions between the scalars and the rest of the fields: 


phi 


— V /1 hi 


^If+r' 1 + —Zf'Z 11 - —tt - —bb 

Vh ^ V h v h v h 


M c _ M r _ 

-CC- TT + 


Vh 


Vh 


a s 

127 TV h 


a 


G%G a ^ + - A UV A^ 


TV h 


pv 1 


+ n i2 hi —- TZ i3 hi Mn " " 1 


v<t, 


NN + V^hihjhk, 


with Vg fc given by 


V-ijk = [^h<t>'TZj2 (VhTZk2 + VpZk l) — Xha^Zj3 {VhTZkZ + ValZkl) 

4 “ 7Zjl ( 6\hVh7Zkl 4 “ ^h<j>V<j>TZk2 ^haVaTZk 3 )] 

+ TZi 2 [Xh^ji ( VhJZk2 + vpZki) — A rpoTZjz (v^TZks + v tT TZk 2 ) 

T 7Zj2 ( 6X ( frV ( pJZk2 4“ Xh(j)Vh^Zkl Xr KJ V(j'JZ-kv ,V 

'TZiZ [ XhalZji {VhJZk3 4” Vq-IZ^y) T X ( - hj 1Z-j 2 {v^lZhS 4“ Vf/TZh^) 
4“ Vj3 ( GX&VfxlZfcs Xho-Vh^Zkl 4“ XfpfjV ( j J lZh2)\ ; 


(3.20) 


(3.21) 


where i,j,k take the values 1,2,3. Note that all scalar vertices containing two or more h 2 's are zero 
due to the Gildener-Weinberg conditions (2.13) and the particular parametrization of the vevs (2.25) 
(see also [58]). Thus, the decay rates for the decays hi —>• h 2 h 2 are zero at tree level. In addition, 
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for all benchmark sets in Table 1, the decay hi —> / 13 A 3 is kinematically forbidden. Therefore, in the 
current framework there are not any lighter scalars that h± can decay to. 

Next, let us consider the total decay widths of all scalars in relation to the corresponding SM Higgs 
total decay width with the same mass T^ ot (Mf l = M^.): 


rtf = nl 


BR^V + BRff + BR™ + BRff 1 + BR|M + BR| m + BR^ 1 + BR§ r M + BR| m 


77 




bb 


xTf 1 (M h = M h 


+ T (hi -» XX) + T(hi -> NN) + T(hi -> hjh k ), (3.22) 

where are the branching ratios of the SM Higgs decays into quarks, leptons or gauge bosons. 

The rest of the decay rates in (3.22) are given by 


T (hi -> XX) = 


327 rvl 


'1 - 


4 Ml 


X 


Ml 


Mi 


K 


1 - 4 7+7 + 12 x 


K 


K. 


I n, 


i2\ 




3/2 


177, 


*3 


r (hi 7 hjh k ) — 


1 - 


2K.+Mn (Ml-Mir 


K. 


+ 


K 


v. 


167T (1 + Sjk) M hi \ 

Using Eq. (3.22), the total decay width Tff of a SM Higgs-like scalar h\ is given as 


ijk 


(3.23) 

(3.24) 

(3.25) 


Tff = cos 2 a cos 2 (3 Tff + T™, (3.26) 

where T^f 1 denotes the total decay width of the SM Higgs with mass My n = 125.09 GeV and T™^ is 
the invisible decay width of the Higgs boson to non-SM states that are kinematically allowed. Namely, 
only when Mx, Mat <62.5 GeV we may have 


r£ v = r (hi -> xx) + r (hi ->• aw) . 


(3.27) 


For completeness, we present in Table 2 the branching ratios of a SM Higgs with mass M ^ = 
125.1 GeV. 


Decay mode 

Branching ratio 

bb 

W+W~ 

99 

T + T ~ 

cc 

Z°Z° 

77 

7 Z° 

0.575 

0.216 

0.0856 

0.0630 

0.0290 

0.0267 

2.28 x 10~ 3 
1.55 x 10 -3 


Table 2: Branching ratios for a SM Higgs boson with M k = 125.1 GeV, for which T® M = 4.08 x 
10 ~ 3 GeV [85]. We did not include the rest of the decay modes because their branching ratios are 
negligible. 
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In order to clarify the deviation of h\ from the SM Higgs, we construct the signal strength parameter 
which can be written as 


= u (jpp hi) BR (hi XX') 
^ hl u SM (pp -s- h) BR sm (h xx) 


(3.28) 


where a, BR are the production cross section and branching ratio of hi and cr SM ,BR SM the corre¬ 
sponding quantities for the SM Higgs. Using (3.26) and a (pp — > h\) = cos 2 a cos 2 (3 u SM {pp —> h), 
the expression (3.28) becomes 

r SM 

M/u = cos 4 a cos 4 P ^5t • (3.29) 

L h 1 

However, due to the smallness of 7£ 2 2 and 7£ 2 3 , the invisible decay width T™ is highly suppressed 
relative to the total decay width and (3.29) simplifies to 


M/u 


~ cos 2 a cos 2 


P- 


(3.30) 


When the Higgs signal strengths from ATLAS and CMS [86-88] are combined [89], one obtains the 
constraint 

p hl > 0.81, @ 95% C.L., (3.31) 

which translates to 

7£n = cos a cos P > 0.9. (3.32) 

Using the benchmark values of the first set of Table 1 we obtain 

Tin = 0.994, (3.33) 


which lies comfortably within the allowed range. Thus, the state h± behaves mostly like the SM 
Higgs boson and is at the moment indistinguishable from it. Run II of the LHC may be able to 
provide a check for the scalar sector of the model if a universal deviation for the SM Higgs couplings 
is established and if new scalar states are discovered. 


4. Dark matter analysis 

Most of the matter and energy content of the Universe still remains a mystery and has therefore been 
dubbed “dark”. Regarding dark matter (DM) in particular, particle physics should be able to explain 
its nature. In the SM alone, however, there is no viable candidate that can play the role of DM. 
Extending the SM, then, becomes again a necessity. On theoretical grounds, the DM particle has 
to be stable and produce the correct relic abundance Ddm- Many direct and indirect DM detection 
experiments are underway or scheduled to become operational in the forthcoming years, aiming also 
at distinguishing between various proposed DM candidates and finally determining its properties. The 
befitting framework for a DM candidate should simultaneously evade limits from previous searches 
and be falsihable in the near future. 

4.1. Boltzmann equation and relic density 

As stated in the Introduction, the rationale behind introducing the hidden or dark SU(2)x gauge sec¬ 
tor is twofold. First, when the new sector is spontaneously broken by means of the Coleman-Weinberg 
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mechanism, the electroweak scale is dynamically generated through the Higgs portal. Second, since 
the SU(2)x gauge symmetry is completely broken by the vev v 0 of the scalar complex doublet 4>, 
the three dark gauge bosons X a acquire equal masses Mx = \gxv<j> and become stable due to a 
remnant global SO(3) symmetry, thus rendering themselves potential WIMP dark matter candidates. 
SU(2)x vector dark matter has been studied in [90-102] and in the context of classical scale invariance 
in [47,48,51,59]. 

After the end of inflation and the assumed subsequent reheating, all particles are in thermal 
equilibrium, while the Universe continues to expand. As the temperature continues to drop, so does 
the interaction rate Tom of the DM particles. Nevertheless, thermal equilibrium cannot be maintained 
and, once Tdm becomes smaller than the expansion rate H of the Universe, a DM number density 
“freeze-out” occurs and dark matter particles decouple from the rest of the light degrees of freedom 
that remain thermalized. Thus, the DM relic abundance survives to the present epoch having the 
value that we observe today. Next, we calculate Ddm, following Refs. [103-105]. 

We start with the Boltzmann equation which describes the evolution of the number density n of 
a given particle species over time. The dark vector bosons X a can both annihilate and semiannihi- 
late [106,107], the relevant processes being listed in Figs. 3-5. 


X Z X W+ X f X N 



X Z X W~ X f X N 


Figure 3: Feynman diagrams for DM annihilation to gauge bosons and fermions. 


X h J X h 5 X h i X hj 



Figure 4: Feynman diagrams for DM annihilation to scalars. 



Figure 5: Feynman diagrams for DM semiannihilation. 


The corresponding Boltzmann equation has the form [59] 


dn 

dt 


+ 3 H n = 


(<Tv) 2 \ 2 ( av )s / \ 

o ( n - n eq )-o n ( n - n eq) 


(4.1) 
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where H is the Hubble expansion parameter, n eq is the number density during equilibrium and (av) 
is the thermally averaged cross section of the DM particles times their relative velocity, with the 
subscripts a and s denoting annihilation and semiannihilation respectively. The thermally averaged 
cross section times velocity is given in the nonrelativistic approximation by [108] 




MJ 


W 


(*) - (2w(s)-u/(s)) 


s=4M| 


(4.2) 


with the quantity w(s) defined as 

™0) = \ i 1 ~ “?) P( s ’ m j’ m k) j d(yC °* 9 ^ ^\M {XX ->■ all) | 2 , (4.3) 

where J0|A1| 2 stands for the matrix element squared of all possible channels, averaging over initial 
polarizations and summing over final spins, f 3 (s,mj,rrik) is the final-state Lorentz invariant phase 
space f 3 (s, rrij, m*,) = ^[l — (mj+mk) 2 /s] 1 ^ 2 [l — (mj—mk) 2 /s ] 1 ^ 2 and s denotes the usual Mandelstam 
variable s = (p\ + P2) 2 = 2 (M| + E1E2 — P1P2 cos 9 ) . Finally, the prime stands for differentiation 
with respect to s/(4M|-) and x is defined as x = Mx/T. The relevant individual cross sections 
necessary for the determination of the total annihilation and semiannihilation cross sections are too 
lengthy to write down. However, these can be found from an analogous calculation in the appendix 
of Ref. [109]. In Fig. 6 we present these cross sections with respect to the dark matter mass M\ for 
a fixed right-handed neutrino mass at Mjy = 240 GeV. 



0 100 200 300 400 500 600 700 800 900 1000 


M x [GeV] 

Figure 6: (color online). This plot shows the thermally averaged total annihilation (purple solid 
line) and semiannihilation (black solid line) cross sections times relative velocity with respect to the 
dark matter mass Mx- The peaks correspond to the poles of the scalar propagators. 


We observe that the thermally averaged semiannihilation cross section is almost an order of mag¬ 
nitude larger than the thermally averaged annihilation cross section. Also, we see two peaks for {crv) a 
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that correspond to Mx = M^/2 and Mx = (for set A in Table 1) and arise due to the form 

of the scalar propagators at s = 4 M\: 


u h, = 


AMl-Ml+iM h V^ 


hi 


hi 


(4.4) 


with T** given in (3.22). There is no peak for the darkon li 2 because its mass varies since it depends 
on Mx- 

Returning to the Boltzmann equation (4.1), it is useful to express it in terms of the comoving 
volume Y = n/s, Y eq = n eq /s, where s is the entropy density, as 


dY 

dx 


-J?T-r4 )-|§(r 2 -ru,), 


3x 2 


s (x 
H(x 


1 ) 


(av) c 


(4.5) 


The entropy density is given by s = 


_ 27Pgr* 
45 


M x 


and the Hubble parameter is given by H = -jf- 


Ml 


45 M P 

in terms of the effective number of relativistic degrees of freedom g * at the time of freeze-out (x = Xf). 
In order to solve (4.5) we may consider the two extreme regions x -C Xf and i » ij, whereupon, 
defining A = Y — Y eq [105] we obtain 


A 


A' 


le? 2 {z a + z s y 

A 2 

— TTT + 2 Z s ) , 


when x <C Xf, 
when x Xf, 


(4.6) 

(4.7) 


where the prime now denotes d/dx. Moreover, if we define A(x/) = cY eq (xf), with c being a constant 
of order one, we can match the solutions of (4.6) and (4.7) and obtain an expression for the freeze-out 
point which can be solved iteratively [105]: 


Xf = In 


0.038 


3 Mx Mp 
yj g*(Tf)xj 


c (c + 2) ( av) a + 2 c (c + 1) (av) 



(4.8) 


We find typical values between Xf ~ 25 — 26 and for the DM mass range that we consider in our 
numerical analysis we use g * = 86.25. Also, for the constant c we use c = 1/2 [110]. The present day 
relic abundance is obtained by integrating (4.7) from x = Xf to x = oo: 


y- 1 = 

OO 


Z a + 2 Z s 

3x 2 


dx. 


(4.9) 


Then, using the mass density of the DM particles today, = MxSocToc and the critical density 
p c = 3H 2 (oo)Mp/(8ir) = 1.054 x 10 -5 h 2 cm' 3 , we hnally obtain the dark matter relic density 


Vt x h 2 


Poo , 2 _o 1-07 X 10 9 GeV- 1 

h — 3 x _ . . ? 

Pc y/fh Mp J (Xf) 



(° v )a + 2 (° v )s 


(4.10) 


where we used = 2891.2 cm 3 for the present day entropy density and h = 0.673 for the Hubble 
scale factor [111]. 

The measured value for the DM relic density is HoMh 2 ± la = 0.1187 ± 0.0017 [111], which is a 
combination of the results from Planck+WP+highL+BAO. In Fig. 7 we scan again the parameter 
space of (gx,Y a ), but this time we also include the points where the DM relic density is saturated 
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within 3 ct (black region). We observe that the DM mass is constrained to be between Mx ~ 710 — 
740 GeV. 


250 



500 1000 1500 2000 2500 


M x [GeV] 

Figure 7: (color online). This plot is the same with Fig. 1 but we also calculated the points where 
the dark gauge bosons can saturate the observed DM relic density at 3cr (black band). 

4.2. Dark matter direct detection 

In recent years, numerous experiments have been set up aiming at directly detecting WIMP dark 
matter. So far these searches have not been fruitful in actually detecting dark matter. However, with 
each new experiment pushing the limits of sensitivity, DM detection could be just around the corner. 

In the present model, the DM candidate X can in principle interact with the nucleons through the 
t-channel exchange of scalar bosons hi . The relevant Feynman diagram is presented in Fig. 8. 



Figure 8 : Feynman diagram for DM-nucleon elastic scattering. 
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This interaction is expressed through the following effective Hamiltonian in the limit of small 
momentum exchange between the DM particle and the nucleon: 


H e g 


V<j> 


X 


i 


'R-rt'Ru 

WT 



(4.11) 


where 77*2 and TZu are the rotation matrix elements from (2.23). The nucleonic matrix element can 
be parametrized as ^ Af \ s ^ jq m q qq\J\f ^ = fM m Mi where mjq- = (m v + m n ) /2 = 0.939 GeV is the 
average nucleon mass and fj\[ = 0.303 [101,112] is the nucleon form factor (see also [113-115]). The 
spin independent dark matter elastic scattering off a nucleon cross section then has the form 


&SI 


Mred 

™h v l 


L\ r M x rn\f ^ 

i 


'R-i2'R-l i 


2 


(4.12) 


where fi re g = Mxmjq-/ (M x + mx) is the DM-nucleon reduced mass, just fi re< i ~ mj\f in our case. 

Employing (4.12), we evaluate the spin independent cross section for various M x and M x masses 
and then, using the experimental results from LUX (2013) [116] and the projected limits from XENON 
IT [117], we construct the plot shown in Fig. 9. 



Figure 9: (color online). The plot shows the DM-nucleon cross section as a function of the DM mass 
for varying M x masses respecting the stability and perturbativity constraints discussed in Sec. 3.1 
(magenta band). The purple solid line corresponds to the experimental limits from LUX (2013) and 
the black solid line corresponds to the anticipated results for XENON IT. 

We find that relatively low M x masses are excluded by LUX (2013). Nevertheless, masses above 
circa 700 GeV, such as those suitable for the saturation of the measured relic density (cf. Fig. 7), 
are favored for detection by XENON IT. 
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5. Summary and Conclusions 


Classical scale symmetry as a framework for model building has recently received a lot of attention, 
mainly due to its appeal as a possible solution to the hierarchy problem through the dynamical 
generation of mass scales. 

In this paper we considered a classically scale invariant version of the Standard Model, enlarged by a 
dark SU(2)x gauge group which incorporates three vector bosons and a scalar field in the fundamental 
representation. We also included a real singlet scalar field and Majorana neutrinos coupled to it. The 
dark sector was radiatively broken through the Coleman-Weinberg mechanism and a mass scale was 
communicated to the electroweak and neutrino sectors through the portal interactions of the dark 
doublet with the Higgs and singlet scalars. 

We started by determining the necessary conditions for the stability of the potential and then 
proceeded in studying the full one-loop scalar potential, employing the Gildener-Weinberg formalism. 
We obtained the scalar masses through a particular parametrization of the scalar vevs, and saw 
that one of these masses, although zero at tree level, received large quantum corrections. Neutrinos 
obtained masses through the realization of a type-I low-energy seesaw mechanism. 

After setting up the model, we proceeded to consider constraints to its set of free parameters 
through stability and perturbativity considerations. Thus, upper limits were obtained for the values 
of the extra gauge coupling gx < 2.51 and the right-handed neutrino Yukawa coupling Y a < 0.31. 
Subsequently, the vevs and the scalar couplings were fixed by the requirement that the correct mass 
for the observed Higgs boson My l] = 125.09 GeV is obtained. Then, the mass of one of the other 
scalars M ^ 3 was readily computed. In order to get the last scalar mass M /, 2 (darkon), we scanned over 
the two-dimensional parameter space of gx = 2Mx/v ( j > and Y a = Mx/v a , since the darkon’s mass 
depends on the masses of every field present in the model (Fig. 1). In addition, we constructed an 
effective Lagrangian describing the interactions of the mixed scalars with the rest of the fields, where 
cubic terms of the form hi — h 2 — /12 were not present due to the Gildener-Weinberg minimization 
conditions (2.13). Since for the benchmark sets of values in Table 1 the mass of /13 turned out to 
be larger than the h\ mass, we concluded that in the current setup the Higgs boson does not decay 
to any lighter states. Afterwards, we calculated the signal strength parameter corresponding to h\ 
and found that the common suppression factor 1Z\i of the Higgs’ couplings with the SM fields is in 
agreement with bounds set by LHC. 

Finally, we identified the extra gauge bosons as WIMP dark matter candidates. We considered 
the Boltzmann equation and solved it semianalytically in the nonrelativistic approximation after 
calculating the total thermally averaged annihilation and semiannihilation cross sections (Fig. 6 ). 
In this way we obtained the dark matter relic density, and by matching it to the observed value we 
constrained the dark matter mass to be in the range Mx ~ 710 — 740 GeV (Fig. 7). Then, considering 
the dark matter elastic scattering off a nucleon we computed the spin-independent scattering cross 
section and compared it with existing and projected limits from direct detection experiments. We 
found that dark matter masses above ~ 700 GeV evade limits set by LUX (2013) but can nevertheless 
be tested in the next years by XENON IT. 

In conclusion, the classically scale invariant model that we considered is a perfectly viable extension 
of the Standard Model, able to dynamically generate the dark matter, neutrino and electroweak scales 
through the multi-Higgs portal while stabilizing the vacuum. It predicts new scalar states that future 
collider searches may be able to discover and predicts vector dark matter with a definite mass range 
that can be probed by direct detection experiments in the years to come. 
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